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Saturn's rings are known to consist of a large number of water 
c-j ice particles. They form a flat disk, as the result of an interplay of 

Q-i angular momentum conservation and the steady loss of energy in 

O dissipative particle collisions. For particles in the size range from a 

few centimeters to about a few meters a power law distribution of 
radii ~ r~ q , with q « 3, is implied by the light scattering properties 
of the rings [1, 2, 3, 4]. In contrast, for larger sizes the distribution 
drops steeply [5, 6, 7, 8] with increasing r. It has been suggested 
^ that this size distribution may arise from a balance between aggre- 

ON gation and fragmentation of ring particles [9, 10, 11, 12, 13, 14, 15] 

but to date neither the power-law dependence, nor the upper size- 
cutoff have been explained or quantified within a unique theory. 
Here we present a new kinetic model for the collisional evolution 
of the size distribution and show that the exponent q is expected 
to be constrained to the interval 2.75 < q < 3.5. An exponential 
cutoff towards larger particle sizes establishes naturally, the cutoff- 
radius being set by the relative frequency of aggregating and dis- 
ruptive collisions. This cutoff is much smaller than the typical scale 
of micro-structure seen in Saturn's rings (~ 100m for self-gravity 
wakes) [16] and our theory represents values averaged over these 
structures. 

A short life-time of Saturn's rings to various evolutionary processes [17, 
18, 19] indicates that ring particles and their size distribution may in fact not 
be primordial but are involved in an active accretion-destruction dynamics [20]. 
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Tidal forces fail to explain the abrupt decay of the size distribution for large 
particles (Supplementary Information) and it is not clear, whether the ob- 
served size distribution is peculiar for Saturn's rings, or if it is universal for 
planetary rings in general. Here we explain and quantify the reported prop- 
erties of the size distribution within a kinetic aggregate model, based solely 
on the hypothesis of binary collisions of ring particles and the classification 
of the possible outcome of collisions as aggregative, restitutive, or disrup- 
tive. We quantify the rate of either collision type, depending on the collision 
speed, using the kinetic theory of granular gases[21] applied to a multi-sort 
system. 

Ring particles may be imagined as "dynamic ephemeral bodies" [11], built 
up of primary grains [13] of a certain size r and mass m . Let the concen- 
tration of ring particles of size k, containing k primary grains, be n^. In 
the kinetic approach (Supplementary Information) these particles are again 
assumed to be spherical, having mass m k = km and radius r k = r /c 1//3 . 
Dense rings composed of hard spheres can be described in the framework of 
the Enskog-Boltzmann theory [22, 23, 24]. In this case the rate of binary colli- 
sions depends on particle dimension and relative velocity. Namely, the cross- 
section for the collision of particles of size i and j scales as r 2 ,^ 1 / 3 + j 1 / 3 ) 2 , 
while the relative speed, being of the order of[4] 0.01 — 0.1 cm/ s, is determined 
by the velocity dispersions (V- 2 ) and (vf) for particles of size % and j, re- 
spectively. The velocity dispersion quantifies the root mean square deviation 
of particle velocities from the orbital speed (~ 20 km/ s). These deviations 
follow a certain distribution, implying a range of inter-particle impact speeds, 
and thus, different collisional outcomes. The detailed analysis of an impact 
(Supplementary Information) shows that for collisions at small relative veloc- 
ity, or small relative kinetic energy E, smaller than a certain threshold energy 
E agg (determined by adhesion alone), particles stick together to form a joint 
aggregate. This occurs because adhesive forces acting between particle sur- 
faces are strong enough to keep them together. For larger velocities particles 
rebound with a partial loss of their kinetic energy. For still larger impact 
speeds the relative kinetic energy exceeds the threshold energy for fragmen- 
tation, -Efrag (again determined by the adhesion only), and particles break 
into pieces. Let be the collision rate (the number of collisions per unit 
time in a unit volume) for particles of size % and j. Only a fraction of these 
collisions, say, A agg z/jj leads to the formation of an aggregate of size (i + j). 
Similarly, only the fraction Af rag z/jj, leads to fragmentation of particles. If the 
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threshold energies and E{ ra , g are some effective constant (Supplementary 
Information), the factors A agg , Af rag and A = Af rag /A agg , characterizing the ra- 
tio of aggregative and disruptive collisions are constant as well. To illustrate 
the approach we consider here a simplified collision model, assuming that 
both particles break completely into monomers in a disruptive collision. A 
more general fragmentation model, with a particular distribution of debris 
size, gives essentially the same result if small fragments dominate (that is, 
size slope of fragments > 4) (Supplementary Information). Schematically, 
we represent restitutive, aggregative and disruptive collisions as 

[«1+b1 H+bl , W+bl [i+j] , [ i] + [i] + ••• + [!]■ 

v v ' 

i+j 

The collision rate ^ oc riifij is proportional to the concentrations of the 
colliding species, rii and rij. Then, with the notations, A agg z/jj = CijUifij and 
Afrag^rj — ^CijUifij, we write the evolution equation following from kinetic 
theory for the concentration of particles of size k > 1, 

dn k 1 00 00 00 

-jr = 7) Y. C ij n i n j ~ C ^ n '^k ~ A Y C ik nin k . (1) 
at 1 i+j=k i=i i=i 

Here the first term in the right hand side describes the rate at which aggre- 
gates of size k are formed in aggregative collisions of particles % and j. By 
conservation of mass, the summation extends over all % and j with i+j — k, 
and the factor | avoids double counting. The other two terms give the rate at 
which the particles of size k disappear in collisions with other particles of any 
size i, due to aggregation (second term) or fragmentation (third term). The 
restitutive collisions do not enter the above equation, since they do not alter 
the concentration of different species. A similar equation may be written for 
the evolution of monomers (Supplementary Information). 

The rate coefficients C^- and the constant A in Eq. (1) are found from 
kinetic theory as functions of particles sizes % and j and threshold energies 
-E agg and -Ef rag (Fig. 1 and Supplementary Information). Moreover, we find 
that these coefficients are self-similar, namely, 

C a i,aj = a ^Cij ) (2) 

which means that re-scaling of particle size by a factor a leaves the form of 
these coefficients unchanged. The constant \i depends on the kinetic model 
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used for the ring particles: If all species have the same average kinetic energy 
(E) = (^TjirLiV^j, the scaling coefficient reads, fi — ^. The opposite case of 
equal velocity dispersion (V 2 ) = (V- 2 ), implies fj, — |. In reality planetary 
rings are in state that is intermediate between these two cases. Indeed, the 
steady state velocity dispersion establishes as a balance of collisional cooling 
and shear heating, leading presumably to a fairly flat decrease of velocity 
dispersions with particle size, rather than to the energy equipartition[25]. 
Hence, ^ < /x < | is expected. 

Systems governed by kinetic equations of the form (1), with self-similar 
coefficients, have the important feature that their solution is mainly deter- 
mined by the degree of homogeneity /i, but not by the particular form [26] 
of Cij. This allows us to use simplified coefficients = C (ij) fl , with the 
same homogeneity exponent /i and to find an analytical solution to Eqs. (I). 
The complete solution is involved (Supplementary Information). However, 
in the limit A <C 1, corresponding physically to rare fragmentation events, 
the steady-state solution has the simple form 

n k ~k-^ 2+ ^e~ x2k /\ (3) 

To obtain the radii distribution f(R), which is constrained by observations, 
we use k ~ r 3 , = R 3 for spherical ring particles and obtain 

f(R)~R-« e -W R ^, g = 5/2 + 3/i, R 3 C = 4r 3 /X 2 . (4) 

For R <C R c the distribution obeys a power-law with the exponent q = 
5/2 + 3/i; for larger radii, R ~ R c , it has an exponential cutoff. 

Two important conclusions follow from this result. First, for the models 
of particle kinetics we study, the power-law exponent q varies within the 
interval 2.75 < q < 3.5. Namely, for the limiting model where the average 
kinetic energy is the same for all particles, independently of their size, fj, — ^, 
and we find q = 2.75. For the model where the velocity dispersion is the same 
for all species, /i — | and q = 3.5. Naturally, these are two extreme cases 
and the ring particles must demonstrate an intermediate behavior, which 
corresponds to the above interval for q. Generally, the size slopes obtained 
from observations [2, 3] consistently lie in this interval, varying between q ohs = 
2.74 and 3. If. Figure I shows a comparison of the theory to the particle 
size distribution derived for a part of Saturn's A ring[2]. Both theoretical 
and observational distributions represent an average over local small-scale 
(~ 100m) ring structure (Supplementary Information). 
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The second conclusion refers to the exponential cutoff of the distribution 
(4). It arises as a direct consequence of the competing fragmentation and 
aggregation processes, so that no other external mechanism is needed. The 
cutoff size, R c is completely determined by the value of A, which depends 
on the aggregation and fragmentation energies -E agg and E{ rag (Fig. 1), and 
the minimal particle size r (Eq. 4). A fit to observations (Fig. 1) gives 
R c = 5.6 m. Using R c , our theory allows to constrain the size r of primary 
grains if the value of the average random velocity V is known (Fig. 1 and 
Supplementary Information). A cutoff radius of R c = 5.6 m is consistent 
with tq in the range of 0.1 — 5 cm for average random velocities V from the 
interval 0.01 — O.lcm/s, again in agreement with the observations[2, 4, 27]. 
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Figure 1: Particle size distribution function. Symbols represent the 
particle size distribution of Saturn's A ring inferred from data obtained by 
the Voyager Radio Science Subsystem (RSS) during a radio occultation of 
the spacecraft by the rings [2]. A fit of the theoretical curve, Eq. (4), is 
shown as a solid line. A cutoff radius of R c = 5.6 m, obtained from the fit, 
corresponds to different combinations of the primary grain size r and average 
random velocity V. Namely, employing Eq. (4) for R c together with A = 

1 — (1 + E a g g / (E))e~ E ^ s ^ E ^ e - E f^ s /( E )^ one obtains a relation between 
ro and V (Supplementary Information). Solving this equation we find that 
r ranges from 4.9 cm to 2 mm, when V varies from 0.01 cm/s to 0.1 cm/s. In 
particular, for V = 0.03 cm/s the primary grain size of r = 1 cm is obtained. 
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1 Collisional mechanisms of particle aggrega- 
tion and fragmentation 

1.1 Derivation of the system of kinetic equations 

In this section we derive the kinetic equations, which describe both aggrega- 
tion and fragmentation of colliding particles. Here we will follow the deriva- 
tion of Ref. [1] done for a similar system, but with another fragmentation 
model. 

Consider a system comprised of particles (monomers) of mass m and 
radius r , which aggregate and form clusters of 2, 3, . . . , k, . . . particles with 
masses = km . For the purpose of kinetic description we assume here 
that all particles are spheres, then the characteristic radius of an agglomerate 
containing k monomers scales with mass as ~ r Q k l / D , where D is the 
dimension of the aggregates. If clusters are compact, D = 3, however, the 
aggregates in principle can be also fractal objects with D < 3. We assume 
that clusters merge when the energy of the relative motion of two colliding 
clusters in the center-of-mass reference frame (the "relative kinetic energy" 
in short) is less than -E agg - In this case a particle of mass (i + j)m is formed. 
If the relative kinetic energy is larger than E agg , but smaller than Ef rag , the 
colliding particles rebound without any change of their properties. Finally, if 
the relative kinetic energy exceeds Ef rag , both colliding particles disintegrate 
into fragments. In what follows we assume that the particles are not too 
large, so that the forces, which keep together particles in aggregates are 
of the adhesive nature, and that the gravitational interactions between the 
constituents may be neglected (see Sec. 1.2.1 for the detailed discussion). 

Note, that recently it has been shown, that all particles with dimensions 
below a certain radius are absent in dense planetary rings [2]. This study 
justify our present model, where no particles with radii r < r are present in 
the system. 

1.1.1 Uniform systems 

To illustrate the main ideas and to introduce the basic notations we first 
consider dilute and spatially uniform systems. Let = f(vi,t) be the mass- 



velocity distribution function which gives the concentration of particles of 
mass rrii with the velocity Vi at time t. The mass-velocity distribution func- 
tion evolves according to the Boltzmann equation 

(v k , t) = /f s + ir + i l r + 4 heat , (i) 

where I k se , I k es and J^ rag are, respectively, the collision integrals describing 
collisions leading to aggregation, restitution, and fragmentation. /£ eat de- 
scribes the viscous heating due to the shearing of the system in its orbital 
motion The first integral reads 

I^(v k ) = \ £ 4 J dv t J dv 3 J deQ (-vtj ■ e) |tiy • e | x (2) 

i+j=k 

x fi (vi) fj (vj) 6 (E agg - E ij )5(m k v k - - rrijVj) 



~ a lj J dv *j J de (-Vkj ■ e ) \v k j ■ e \ 
j 



x 



x fk (v k ) fj (vj) 6 (E agg - E kj ) . 

Here = r (i 1 / 3 + j 1 / 3 ) is the sum of radii of the two clusters, while 
Tfik = ™>i + rrij and m k v k = rriiVi + rrijVj, due to the conservation of mass 
and momentum. The unit vector e specifies the direction of the inter-center 
vector at the collision instant. We also introduce the reduced 
mimj/(mi + rrij), the relative velocity, Vij = Vi — Vj, and the kinetic energy 
of the relative motion, = \\iijv\j. The factors in the integrand in Eq. 2 
have their usual meaning (see e.g. [3]): afj ■ e\ defines the volume of 
the collision cylinder, (—Vij ■ e ) selects only approaching particles, while 
(-Eagg — Eij) guarantees that the relative kinetic energy does not exceed 
-E agg to cause the aggregation. The first sum in the right-hand side of Eq. 2 
refers to collisions where a cluster of mass k is formed from smaller clusters 
of masses i and j, while the second sum describes the collisions of /c-clusters 
with all other aggregates. 

For collisions leading to fragmentation we have 

4 rag K) = g E °ij J dv *j J dv *i J deQ {-Vij ■ e) |{% • e\ fj (vj) f\ (^) x 

i j 

x 6 (E^ - £ frag ) (q ki (v k , Vi, Vj) + q kj (v k , v { , vj)) (3) 
- (! - $k,i) J d^ J deQ (-v ki ■ e) \v ki ■ e \ x 
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where we define the kinetic energy of the relative normal motion, Ef- = 
\lHj{vij • e) 2 , with {vij • e) being the normal relative velocity. Note that in 
contrast to the case of aggregation, where both normal and tangential com- 
ponents must be small, so that the total energy of the relative motion E^ 
matters, for a fragmentation process only the relative normal motion is im- 
portant: It causes a compression and the subsequent breakage of particles' 
material. Hence the kinetic energy of the relative normal motion, Efj must 
exceed some threshold. The first sum in Eq. 3 describes the collision of par- 
ticles i and j with the relative kinetic energy of the normal motion above the 
fragmentation threshold Ef rag and qki{vk, Vi, Vj) gives the number of debris 
of mass mfc = mo k with the velocity Vk, when a particle of mass rrii = m i 
desintegrates at a collision with a particle of mass m,j = moj, provided that 
the according pre-collision velocities are Vi and Vj. Obviously, = if 
k > %. The function qki{vk, fy, Vj) depends on the particular collision model. 
For the purpose of the present study we restrict the analytical model on the 
case of complete disruption of clusters into monomers (section 2). However, 
numerical results for more general cases, when particles decay in debris with 
a power-law size distribution, are also shown for comparison in section 2.2. 
The second sum describes the process, when particles of mass k ^ 1 break 
in collisions with all other particles. The collision integral for the restitutive 
collisions has its usual form [3], 

ir(v k ) = £ °li J d % J • e ") Wki ■ e| x (4) 

i 

x (e- kt 2 fk K') fi «') - fk (vk) fi («*)) © (Eh - E agg ) 9 (E irag - E ki ) , 

with a slight modification to account the requirements for the relative kinetic 
energy: Eij > E agg and Ef- < E {rag . Here vj! and v" are the initial velocities 
for the so-called inverse collision, which results with the velocities Vk and v*i, 
and Eki is the restitution coefficient for a collision between particles of size k 
and i, e.g. [3]. 

For the heating term we adopt the following simple model, e.g. [4] 

4 heat = y|U(^), (5) 
where T k is the rate of heating of the species with mass m k . 
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Owing to the permanent aggregation and fragmentation we have a dy- 
namically sustained mixture of particles of different mass. Mixtures of dissi- 
pative particles, generally, have different velocity dispersion, or mean kinetic 
energy ("granular temperature") of each species. We define the partial num- 
ber density rii and mean kinetic energy (Ei) of the species as [5], 

rn = J dvifi(vi) , ^(Ei) = J ^m^ 2 f{m u Vi)dvi , (6) 

We assume that the distribution function fi(vi,i) may be written as [6, 7, 5] 

fi(i>i, t) = H ^f 4>i{ci) Ci = — , (7) 

v 0,i\ T ) V 0,i 

where Vq j(t) = 2{Ej){t) /rrii is the thermal velocity and 0(q) the reduced dis- 
tribution function. For the force-free granular mixtures the velocity distribu- 
tion functions of the components are not far from the Maxwellian distribution 
[5], which reads in terms of the reduced velocity c — v/v , 

<M3=7r- 3 / 2 exp(-c 2 ). (8) 

The equipartition between different components, however, breaks down in 
general for granular systems, in the sense that the average kinetic energies 
for each size group (Ei) are not equal [5]. 

Integrating Eq. 1 over v k we obtain the equations for the zero-order mo- 
ments of the velocity distribution functions fk, that is, for the concentra- 
tions Uk- Taking into account that collisions resulting in rebounds do not 
change the concentrations of different species and using Eqs. 2-3 together 
with Eqs. 6-8 we arrive at the rate equations 

^-n k = \ C ij n i n j ~ n kJ2 C ki n i - A kiKkni (1 - 6 kl ) 

dt 1 i+j=k i i 

+ \ A ij n i n j ( x k(i) + x k (j)) , (9) 

where x k (i) is the total number of debris of size k, produced in the disruption 
of a particle of size i. 1 Note that x k (i) — if k > i and J2kk x k(i) = i- In 



1 The value of Xk{i) may be expressed in terms of 
Qk,i,j{z) — Qk,i,j((vij ■ e) 2 ) — J qki(vk,v i} Vj)dvk, where we assume that the disruption 
outcome depends only on the square of relative normal velocity z — (vij ■ e} 2 as Xk(i) = 
Xk(i,j) - luijBije^^ / °° e-^ B »'Q k ,i,j{z)Q(z -2E bae /m j )dz . 
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what follows we will use the model expressions for the coefficients Xk(i), like 
x k(i) — ifiik, or x k{i) ~ k~ a , which do not depend on the velocity dispersion 
or aggregation and fragmentation energies. 

The rate coefficients Cy and A^ in the above equation read: 



Cij 



v ij Aagg 



Vij A 



frag 



V, 



l 3 



24 



\ 




(10) 



A agg = 1 - (1 + BijEzgg) exp (-BijEggg) 



Afrag = exp 1 



-BijEft 



5 ) 



B 



m,i + rrij 
(Ei)rrij + (Ej)mi ' 



Here VijUiiij give the collision frequency between particles of size i and j 
and A agg and Af rag give respectively the fraction of aggregative and disruptive 
collisions. It is useful to verify that the above kinetic equation (9) fulfills the 
condition of mass conservation, Y.k^ m o n k — M — const., where M is the 
total mass density. 

If -Ef m g = const., -E ag g = const, and (Ej) = (E) = E av , the rate coefficients 
A^ and are proportional and differ only by the multiplicative constant 

A Af r ag/ ^agg- 



2at\ 2ttE r 



Aij XCij , 



(rrij 1 + rrij x ) 



1 — (1 + E &gz /E & 



g--Eagg 



A = 



1 - (1 + £ ag g/£av) e~ E ^ 



1 



In what follows we assume for simplicity that the above form of the kinetic 
coefficients holds generally, that is, the dependence of the rate kernel on the 
particles size occurs only through the collision frequency 



a 



const, x v. 



VI 



AC, 



13 1 



(12) 



where A is a constant. 



1.1.2 Self-gravity wakes and ring structure 

Saturn's rings are not uniform but exhibit a large variety of structures [8]. 
One example are the self-gravity wakes [9, 10, 11, 12], arising from self- 
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gravitational instability, forming a transient and fluctuating pattern in the 
surface mass density of the rings, canted relative to the azimuthal direc- 
tion, with a typical length scale of about L w ~ 10 2 m, one Toomre critical 
wavenlength [13]. Thus, to describe adequately particle kinetics one needs 
in principle to take into account effects of dense packing, non-homogeneity, 
and gravitational interactions between particles. Two important comments 
are to be done in this respect. First, due to the low velocity dispersion of 
particles in the dense parts of the wakes, the collision duration is still sig- 
nificantly smaller than the time between particle collisions. This implies the 
validity of the assumption of binary collisions, as the dominant mechanism 
of particles' kinetics. For the same reason, the the method of molecular dy- 
namics simulations, based on the binary collision hypothesis, is applicable 
to study the dynamics of dense rings [14]. Therefore a kinetic description 
in terms of the Enskog-Boltzmann equation is possible [15]. Although this 
Markovian equation ignores memory effects in particle kinetics, it may be 
still applicable, when the mean free path is comparable to, or even smaller 
then the particle size [16]. Second, the characteristic length scale of the den- 
sity wakes, L w , and the upper cut-off radius (about 10m) are well separated. 
This allows to neglect variations of density and distribution functions on the 
latter length scale and use a local approximation for the distribution function 
of two particles at contact: 

f2(vk,r- er k ,vi,f+ eri,t) ~ g2(cri k )f k (v k ,r,t)fi(v l ,r,t) (13) 

where f'2 is the two-particle distribution function, for particles of radius r k , 
and ri, which have a contact at point r, the unit vector e joins the centers 
of particles and (fe^/fc) is the contact value of the pair distribution function, 
which may be well approximated by its equilibrium value for the hard sphere 
fluid; explicit expressions for g2{&i k ) m ay be found, e.g. in [3]. In the local 
approximation (Eq. 13) the collision integrals depend only on local values at 
a particular space point r. This significantly simplifies the kinetic description 
of a high-density gas, since the density effects are taken into account in this 
approach by the multiplicative Enskog factor, #2(0"^), leaving the structure of 
the collision integrals unchanged. Therefore the Enskog-Boltzmann equation 
valid for the case of wakes reads: 

d - d 

-Q t fk (v k , r, t)+v k - V/fc (v k , f, t) + F k (f) ■ — f k (v k , f, t) = (14) 

= /f s (f) + /r(f) + /i ras (f), 
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where f k (v k ,f,t) is the velocity distribution function of particles of size k 
(comprised of k- monomers) , which depends on the space coordinate r and 

— * 

F k {r) is the total gravitational force, acting on the particle of size k, which 
includes both the gravitational force from the central planet as well as self- 
gravitation of the ring particles. In what follows we do not need an explicit 
expression for this term. The collision integrals in the r.h.s. of Eq. 14 have 
the same form as in the previous case of a uniform system, with the only 
difference, that they depend on local parameters taken at a point r, and that 
the collision cross-sections are re-normalized according to the rule: 

°ij — >■ ^ij^ij) (15) 

which accounts for the high-density effects in the local approximation, see 
e.g. [3]. Note that we do not need here the term I k eat (r), which mimicked 
heating for the model of a uniform gas, since the Eq. 14 implicitly contains 
the spacial gradients and fluxes, responsible for the heating. 

If we now integrate (as in the case of a uniform gas) the kinetic equation 14 
over v k , we find the equations for the concentrations n k (r): 

d - - 1 

— n k (r) + V • jfc(f) = - ^ Cijirjniirjrijir) - n k (f) ^ C ki (r)ni(r) (16) 

i+j=k i 

+ ^A ki (f)n k {f)n i {f) (1 - S kl ) + ^^Aijif)ni{f)nj{r) (x k (i) + x k (j)) , 

i i,j 

where the partial flux jkif) is defined as 

Jk(r) = J v k f k (vk, r, t) dv k . 

It describes the macroscopic (hydrodynamic flux), associated with particles 
of size k. The important difference of Eq. 16 from Eq. 9, used for a uniform 
system, is the spatial dependence of the kinetic coefficients Aij and CV,-. 
Although the structure of these coefficients coincides with that of Eqs.10 
(apart from the trivial re-normalization, Eq. 15), all quantities here are local. 
Naturally, the local velocity dispersion [E^if) in the dense parts of the wakes 
significantly differs from that of the dilute regions inbetween. 

Now we average Eqs. 16 over a suitable control volume V, which contains 
a large number of wakes. Applying then Green's theorem and taking into 
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account that the surface integral of the flux j k over the boundary of the 
control volume vanishes as 1/V, 



J V • j k df = J j k ■ ds = 0, 



we finally arrive at a set of equations of the same form as Eqs. 9, but with 
the space-averaged quantities: 

^rik = 7, Cijfiifij -n k J2 CkiUi - A ki n k ni (1 - S kl ) 

i+j=k i i 

+ IY, AijUiUj (x k (i) + x k (j)) . (17) 
Here, by the definition, 

n k = — J n k (f)df. 

and 

1 1 f 

Clj = f^n v I n i(^ n j(^ C ij(^ d? ^ 

where Cij(f) are defined by Eqs. 10 with the local velocity dispersions (Ei) (r). 
Similar expression holds true for the coefficients Aij. 

It is important to note that the coefficients and are density- 
weighted quantities. Therefore the contribution to the average value is pro- 
portional to the local density. This in turn implies that the values of these 
coefficients practically coincide with these for the dense part of the wakes, 

C< _ /~i{dense part) T t {dense part) 

— ^ij Aij — 

Hence we conclude that the kinetic equations for the average concentra- 
tions of particles n k coincide with the previously derived equations for n k for 
the case of a uniform system. In what follows we will use n k , Cij and Ay for 
the notation brevity, keeping although in mind that they correspond to the 
average values, that are almost equal to these values in the dense part of the 
wakes. 
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1.2 Tidal fores and gravitational versus adhesive en- 
ergy of aggregates 

1.2.1 Impact of the tidal forces on the aggregate disruption 

In the previous section we considered only the collisional mechanism of ag- 
gregates disruption. The rings particles, however, also suffer from tidal forces 
and we wish to estimate the impact of these forces on the aggregates stabil- 
ity. If the critical stress of the material of a body subjected to tidal forces is 
P max , the condition of the tidal disruption reads [17], 

P max < 1.68 p 6 P^ 2 (18) 

where pb is the density of the body, Rb its radius and Q is the Kepler fre- 
quency. For solid ice the above equation yields 200 km for the critical radius 
of tidal disruption. In application to aggregates of radius R comprised of 
spherical ice particles of radius r , we can estimate the critical pressure and 
critical radius as follows. The area of the equatorial aggregate cross-section 
irR 2 contains 4>nR 2 /-nr^ = 0(P/r o ) 2 particles' contacts, where is the mean 
packing fraction of the aggregate composed of primary grains. All primary 
particles are kept in the aggregate by the adhesive forces. According to the 
JKR theory [18] the maximal pulling force which an adhesive contact can 
resist, reads, 

fo = \m r -% (19) 

where 7 is the surface tension of the particle material. The total force which 
can resist a body with the cross-section nR 2 is (p(R/r Q ) 2 f , hence the maximal 
resistance to pulling may be estimated as 

Pmax = (20) 

4 r 

Using r = 1cm, 7 = 0.74 N/m and <p = 0.4, we obtain P max = 22.2N/m 2 . 
With the density of the aggregates of p ice = 400 kg/m 3 and VI = 1.5-10~ 4 s _1 
we obtain the disruption radius of 1.211 • 10 3 m = 1.2 km. If we assume 
that the primary grain size is ten times bigger, ro = 10 cm, and the surface 
tension is 10 times smaller, due to the surface roughness of the particles, 
7 = 0.074 N/m, we still obtain the disruption radius of 121m, which is one 
order of magnitude larger than the cutoff radius R c , observed for the particle 
size distribution in the Rings. Hence we conclude that the disruption of 
particles due to the tidal forces may be safely neglected. 
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1.2.2 Gravitational versus adhesive energy of aggregates 

In our study we concentrate on the aggregates, comprised of particles kept by 
adhesive forces. It is instructive to compare the total energy of all adhesive 
bonds with the gravitational energy of aggregates (for detailed discussion of 
this issue see the forthcoming article [19]). The total gravitational energy 
reads, 

Gi, mi m 3 G r [ „ p{?i)p(?2) 3GM 2 



2 f-f . f-f \fi — fj\ 2 J J \f\ — r*2 5R 

i=i j=i,j^i \ 1 j\ ' ' 

where G is the gravitational constant and we approximate the discrete dis- 
tribution of iV particles of mass rrii at spatial points r*j by the continuum 
uniform mass distribution with density p{r) = Pi ce (f>; here p icc ~ 10 3 kg/m 3 
is the density of the constituents (ice particles) and <fi ~ 0.4 is the average 
packing fraction of the aggregate. The total mass of the aggregate of radius 
R then reads, 

N 4 

M = J2rrii = -7cR 3 p icc (j) . 
i=i 6 

This gives for the gravitational energy in the SI units: 

Eq = 1.124 - i? 5 [J] . 

The total adhesive energy, equal to the energy of all adhesive bonds in ag- 
gregate of radius R, comprised of spheres of radius r , is computed below 
in section 4. The result reads in SI units (this energy is termed as "average 
fragmentation energy" below) 

S adh = 8.2.10- 7 i2 3 ro 5/3 [J] 

For the typical values of R ~ lm and 2tq ~ 0.01m, the adhesive energy 
is almost two order of magnitude larger than the gravitational energy. For 
larger radii R, close to the cut-off values, R 5m, the adhesive energy is 
still few times larger than the gravitational one. 

These estimates motivate our approximate treatment of the aggregates: 
With an acceptable accuracy, we will neglect the gravitational interactions 
of particles in the aggregates and take into account only the adhesive inter- 
actions. 
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2 Solution of the rate equations 



2.1 Analytical approach 

We start the analysis from the simplest fragmentation model - the complete 
fragmentation into monomers of both colliding particles; in this case Xk(i) = 
i5ik- We also assume that the mean kinetic energy is the same for all species, 
so that Aij = XCij. In this case the rate equations take the simplified form, 

dr rr = \ dj n i n j - (1 + A W C kjnj (21) 

al Z i+j=k j>l 

for evolution of clusters of mass k > 1. The parameter A, given by equa- 
tion (11), characterizes the relative weight of fragmentation. The governing 
equation for the monomer density reads 

^jf = ~ni Yl C ^ n 3 + ^ E °ij (* + j) n i n i + An i J C ij n j ( 22 ) 

al j>l Z i,j>2 j>2 

Similar equations with terms responsible for the spontaneous fragmentation 
of aggregates have been analyzed in Ref. [20, 21] in the context of rain drop 
formation. 



2.1.1 Constant kernel 

Evolution kinetics. Consider first a model with mass-independent kernel 
Cij = C . For this kernel, Eqs. 21-22 simplify to 

= \ E run, - (1 + \)n k N (23) 

i+j=k 

and 

^1 = - ni N + A(l - m)N . (24) 

where we re-scale time t — > Con t and densities n& — > n^n^, leaving for the 
re-scaled quantities the old notations. Here 

N(t) = E n M (25) 

is the total cluster density. n$ in the last equations (23)- (25) is the character- 
istic concentration of particles, say the concentration of monomers, and C is 
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the characteristic value of the (equal) rate coefficients CV,-, so that CqUq gives 
the characteristic frequency of aggregative collisions. Summing all equations 
(23) and equation (24) we arrive at a closed equation for the cluster density: 

^ = -N 2 + 2A(1 - N)N (26) 
at 

Consider for concreteness the mono-disperse initial condition n^(t = 0) = 
5k,i- Solving (26) subject to N(0) = 1 gives 

JV W = i + 2A-e-« - <27) 

Plugging this into (24) and solving subject to ni(0) = 1 we obtain n\{t): 



ni(t) = ^ + _L_ A (t) with A(t) 



(1 + 2A)e 2A< - 1 
2A 



2(1+A) 
1+2A 



(28) 

and then, recursively, from equation (23) all rik(t) with k > 1. These rela- 
tions, however, quickly get very unwieldy as the mass k grows and we do 
not present them here. The straightforward analysis shows that the system 
reaches a non-trivial steady state, n^(t) — > rif, and the characteristic relax- 
ation time is essentially the relaxation time for the total density of clusters 
N(t). From Eq. 27 follows the dimensionless relaxation time, t~1 = 2 A; 
therefore, from Eq. 11 we obtain the dimensional relaxation time 

r-, 1 ~ C n A ~ RiVn e- E ^ E ™ , 

which actually coincides with the average collision time for disruptive col- 
lisions. Here i? av is the characteristic size of particles and V is their mean 
random velocity. With the mean packing fraction u = (47r/3)-R^ v n for the 
ring particles, the last equation takes the form, 

r^ 1 ~iWC 1 e~ Wi5av (29) 

Using V = 1CT 4 — 10 _3 m/s, the average particle size _R av ~lm, the packing 
fraction u = 1CT 2 — 10 _1 and the exponential factor, estimated for this 
random velocity interval, as 4 x 1CT 14 — 3 x 1CT 7 (see Section 4), we obtain 
as a lower boundary for relaxation time of the system, 

r re i > 1, 000 years . 
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This is the characteristic time for the rings to return from a perturbed state 
to their proper steady-state size distribution. 

Steady- state size distribution. To obtain the steady-state solution we 
write the governing equations 



0= ntrij - (1 + X)n k N 

i+j=k 



(30) 



which suggests to use the generating function technique. Multiplying (30) 
by z k and summing over k — 2, 3, . . . we obtain a simple quadratic equation 



M{zf - 2(1 + X)NAf(z) + 2(1 + \)Nmz = 
for the generating function 

N{z) = J2n k z k 

k>l 

Solving (31) yields 

N{z) = (1 + A)7V 

Expanding Af(z) we arrive at 

N 



(31) 



(32) 



l-Wi- 



2ni 



;i + a)a^ 



(33) 



n k = 



;i + a) 



2ni 



1 k 



(1 + X)N 



r(fc-|) 

r(A; + l) 



(34) 



Using the already known expressions (27, 28) for the (steady state) densities 
of the monomers and the total number of clusters 



A 



2A 



ni= lTA' "=1^ 



(35) 



we recast (34) into the form 



n k = Aa kT{k $ 



where 



1 A(l + A) 
1 + 2A 



T(A; + 1) 



, a(X) = 1 - 



1 + A^ 



(36) 



(37) 
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When Acl, equation (36) becomes 



nk ~7^ rx^TT) (38) 



which can be further simplified to 



rik = e 

'7T 



A -* k k-V 2 (39) 



when k ^> 1. Thus in the region k < A 2 , the mass distribution exhibits a 
power-law dependence with an exponential cutoff for larger k. 



2.1.2 Generalized Product Kernels 

Equal mean energy. We again assume that Ef r3ug = const., E agg = const, 
and that all species have the same mean energy E av . Then the rate constants 
read, 

Cij = c (i 1/3 + j 1/3 f (r 1 + r 1 ) V2 > = XC v ( 4 °) 

where the constants Cq and A are determined by Eqs. 11. 

It seems impossible to solve (21) for the rate coefficients (40); even in the 
steady state, the equations appear intractable. An old trick in the field of 
aggregation is to replace the kernel by a simpler one with the same degree 
of homogeneity. The reason is that the homogeneity degree is often the 
only characteristics that affects the behavior. Therefore one can use the 
generalized product kernel 

ICa = (ijY (41) 
with the same degree of the homogeneity 

2/1 = 2/3-1/2 = 1/6, or fx = 1/12 

which is analytically solvable. Indeed, the governing equations in this case 
become (we again use the re-scaled densities nu — > nonk), 

= E - (! + x ) l *M (42) 

i+j=k 

in the steady state. Here we have used the notation 

l k = k»n k , M = E4 (43) 
k>i 
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Equation (42) are mathematically identical to Eqs. 30 and therefore they 
admit the same solution (recall that l\ = rii) 



M{z) = (1 + A)M 



1 



2m 



(1 + \)M 



where JA(z) is the generating function 



k>l 



(44) 



(45) 



Therefore the solution reads 

M 



47T 



1 + A) 



2m 
_(1 + A)M 



r(fc-|) 

k»-F(k + l) 



(46) 



The analysis shows that 2ni/[(l + A)M] = a with a given by (37) does not 
depend on \i. Thus the solution (46) can be re-written as 



;i + a)m , r(fc-|) 



2v^r k»-r(k + l) 



(47) 



The quantity M can in principle be found from the sum rule and the behavior 
again simplifies in the small A limit. Asymptotically (when 1 <C k < A~ 2 ) 
we get 



M 



(48) 



which is again the power-law behavior with the exponential cutoff. Note that 
the scaling, n k ~ fc -3 / 2- ^ has been obtained previously in Ref. [21]. 

Equal mean dispersion velocity. If we now assume that all ring particles 
have the same mean velocity dispersion, (V 2 ) = (V 2 ), that is (Ei) = rrii (V 2 ), 
Eqs. 10 give 



C ij = C (iW+j 1 ' 3 ) 2 , 
with the homogeneity degree 
2/x = 2/3 , 



A 



XCij , 



(49) 



or 



A* = V3, 



different from that obtained for the model of equal mean kinetic energy of 
rings particles. Note that the model of equal velocity dispersion implies 
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that (Ei) ~ i, while for the model of equal mean energy [Ej) ~ i°. For 
the former case fi = 1/12, and for the latter /i = 1/3. Physically, one 
expects that the energy of ring particles obeys the relation, (E^ ~ i 1 , where 
< 7 < 1. Indeed, 7 can not be negative, since this would imply that very big 
particles do not have velocity dispersion, which is only possible for unrealistic 
condition of the collision-free motion. At the same time the condition 7 > 1 
is also not realistic, since it would imply that the energy, pumped in due to 
collisions, increases with particle mass faster than linearly; obviously, there is 
no evidence of such processes. Hence we conclude that 7 is limited within the 
interval [0, 1], and, correspondingly, \i obeys the restrictions 1/12</i<1/3. 

2.2 Numerical analysis of the rate equations 

To verify the validity of various approximations made in the analytical so- 
lution resented in the previous section we additionally perform a numerical 
analysis. We explicitly solve the system of rate equations. First, we observe 
that the exact rate coefficients and the rate coefficients for the product kernel 
with the same degree of homogeneity indeed give practically identical results, 
(Fig. 1). We also compare the fragmentation model with complete breakage 
(studied in the previous section) to a model with power law distributed de- 
bris sizes (xk(i) ~ k~ a ). The results are shown in Fig. 2. As one can see from 
the figure the both models (of the complete disruption and of the breakage 
with the power-law distribution of debris size) develop practically the same 
slope if the distribution of debris size in a collision is steep enough, a > 2. 
The steep distribution implies strong domination of small fragments, which 
seems very plausible from the physical point of view: Indeed, the aggregates 
are relatively loose objects, with the small average coordination number. For 
such object it would be much easier to create small debris in a fragmentation 
process, as compared to the case of large coordination number. 

3 Size Distribution 

In observations the size distribution of particles can be constrained. Using the 
relation k ~ R 3 (for spherical particles) in conjunction with n^dk = f(R)dR 
we find that n& ~ k~ 3 / 2 ~^ implies 

f(R) ~ R-Ve-W*** 3 g = 5/2 + 3/i R 3 C = 4r 3 /A 2 (50) 
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k 



Figure 1: Size distribution from the model given by equations (21, 22) for 
two different kernels with the same degree of homogeneity - the ballistic 
kernel and the simplified product kernel. The slope of the power law regime 
of the theoretical prediction, n& oc k 19 ^ 12 , is shown as a dotted line. The 
fragmentation parameter is A = 0.02. 
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Figure 2: Size distributions for models with a power law distribution of 
debris sizes ~ k~ a with a = 10, 5, 3, 2, and 1.5 for the case of constant 
reaction kernel, /i = 0. For steep size distributions of the debris, a > 2, the 
resultant size distribution of aggregates shows practically the same power law 
and exponential cutoff, nu ~ k~ 3 ^ 2 exp(— x k), as the model with a complete 
disruption into monomers. The parameter x, obtained by fitting, coincides 
for all a = 10, 5, 3, and starts to deviate for the frontier value a = 2, 
where the exponent of the power-law also starts to differ from 3/2, taking 
successively descending values. 
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For the physically relevant distribution in three dimensions the appropriate 
generalized product kernel (41) has fi = 1/12 and Eq. 50 implies f(R) ~ 
.R-ii/4. This prediction is in excellent agreement with measurements [22]. 

As it has been shown above, the exponent fi can vary within the interval, 
1/12 < fi < 1/3, which implies that the exponent q for the radii distribution 
obeys the condition, 2.75 < q < 3.5. Occultation experiments give values for 
q in the range from 2.71 to 3.11 [22, 23]. 



4 Calculation of aggregation and fragmenta- 
tion energies 



Here we present estimates for E &gg and -Ef rag . We assume that these quanti- 
ties are some effective constants which allow an adequate description of the 
system kinetics in average. All constants below are expressed in SI units and 
for the notation brevity we skip their dimension. 

To estimate the aggregation energy we use the result of [24] for collisions 
of particles with adhesion. The threshold energy of aggregation reads in this 
case: 



Here go — 1-457 is a constant, 7 is the surface tension, D = (3/2) (1 — u 2 )/Y, 
where Y and v are respectively the Young modulus and Poisson ratio and 
R e g = R 1 R 2 /(Ri + #2)- Using material parameters for ice, Y = 7 x 10 9 Pa, 
v = 0.25, 7 = 0.74 N/m, we obtain E &gg = 8.72 x lCT 7 /?^ 3 J, where i? av is 
the mean particles radius. The corresponding mean kinetic energy reads, 



Here i? av is the mean aggregate size, p = 10 3 kg/m 3 is the ice density and 
gives the mean packing fraction of particles, so that (f>p is the actual density 
of the aggregate. With ~ 0.4 we obtain, E av = 837.8 Rl v V 2 J. 

The above Eq. 51 indicates the energy required to separate particles 
bounded by one adhesive contact. To break a contact of two monomers 
of radius r the energy of E = E agg (R eS = r /2) is needed. To destroy com- 
pletely an aggregate which has a total number of N c adhesive contacts, one 
needs the energy Ef va , g = N c Eq. The number of contacts equals the number 
of monomers in the aggregate of size R, given by (4/3)7ri? 3 0/((4/3)7rrQ) = 




(51) 




(52) 
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0(i?/r o ) 3 , times half the average coordination number n c . As a result we 
obtain for the average fragmentation energy, 

E bag = 2' 7 / 3 q (tt^D 2 ) 1/3 0n c (R av /r o ) 3 . (53) 

This yields for n c = 4.7, found for the random packing of spheres [25, 26] 
E bag = 8.2 x 10- 7 R 3 av ro 5/3 J. Finally, using Eq. 11, we obtain for the coeffi- 
cient A: 



A 



1-1 + 



1.04 x lO -9 ' 



V 2 R, 



5/3 



exp 



9.8 x 10~ 10 



V 2 r° 



5/3 



(54) 



The value of i? av is to be obtained self-consistently. Namely, we use the 
so-called " mass- averaged" radius 



Ra 



J™f(R)R±dR 
J r °°f(R)R*dR' 



which for the distribution (50) yields, 



R„ 



r(*?)J 



R c ) 



(55) 



where we neglect small terms of the order of (r /i? c ) 4 . For q 3 we obtain, 
R av ~ [r(|)/r(|)l R c ~ 0.506 R c . Using Eq. 50, which relates r and A, 
the value of R c = 5.5 m for the cutoff radius, obtained by fitting to the 
observational data, and the mean "mass-averaged" radius R av = 0.506 R c = 
2.79 m, we obtain the following equation, which constrains the monomer 
radius r for different values of the average random velocity V: 



r 3 ^ 2 = A exp {aV 2 r 5 ^ 3 ) 



A = 



- A exp [av r j 
1 - (l + bV- 2 ) exp(-bV- 2 ) 



(56) 



where a = 9.78 x 10" 10 , b = 1.88 x 10~ 10 and c = 0.155. Hence Eq. 56 
determines the radius of a primary grain for each value of average random 
velocity V. Solving this equation we find that r ranges from 4.9 cm to 2 mm, 
when the random velocity V varies from 0.01 cm/s to 0.1 cm/s. In particular, 
the random velocity V = 0.03 cm/s implies a primary grain size of r = 1 cm. 
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Figure 3: The dependence of the primary grain size tq on the chosen value 
of the average random velocity V. The calculations have been done using 
Eqs. 56 for R c = 5.5m (solid line), R c = 2m (squares) and R c = 15m 
(circles). As it may be seen from the figure, different cutoff radii correspond 
to the same combinations of possible primary grain size and average random 
velocity. Note that R c and V are in fact related additionally by the energy 
budget of the ring. 

Similar calculations may be done for the other value of R c - in this case 
the constants in Eq. 56 read, a = 9.78 x 1(T 10 , b = 3.23 x 1(T 9 R~ 5/3 and 
c = 2R~ 3 / 2 . 

Fig. 3 demonstrates the dependence of the primary grain size r on the 
chosen value of the average random velocity V. The calculations have been 
done for R c = 5.5 m, obtained by fitting to the observational data of [22] and 
for the two other values of the cutoff radius, R c = 2m and R c = 15 m. It may 
be seen from the figure that the cutoff radius has practically no impact on 
the value of r for a given velocity V. Now we can estimate the exponential 
factor in Eq. 29 of Sec. 2.1.1. Using Eqs. 52 and (53) we find 



where a = 9.78 x 1CT 10 and ro is to be taken as a function of V from the 
equation (56). Hence this factor strongly depends on the chosen value of 
random velocity. For the random velocity interval of 10™ 4 — 10~ 3 m/s we 



exp (-£ frag /£ av ) = exp ( 



—ar. 



o 
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obtain, 



4.2 x 10" 14 < exp 



a (r (y)) 5/3 V- 2 < 3.4 x 10~ 7 , 



that is, 



4.2 x 1(T 14 < exp (-E iTag /E av ) < 3.4 x 1(T 7 . 



5 The dependence of the upper cutoff radius 
R c on the distance from the planet 

Now we discuss the dependence of the upper cutoff radius R c on the distance 
from Saturn r. According to Eqs. 50 and f f 



R r 



22/3^-2/3 = 2 2/3 



r 



1 - (1 + E agg /E : 



2/3 



,2-Bfrag/S-Ba 



Using E agg and E av found in the previous Section 4 we obtain, 
E agg /E av = 3.26 x 1(T 9 V~ 2 R: 5/3 , 



(57) 



(58) 



which shows that for 1(T 4 < V < 10- 3 m/s [27, 28] and 2.4 < R c < 11.2 m 
[22] the ratio E^/E^ is constrained to the interval 0.006 < E agg / E av < 
0.075, that is this ratio is small. Hence one can use in Eq. 57 an expansion 
in E agg /E w which yields, 



Rc = 2 2 / 3 r (£ agg /£ av ) 4/3 e 2 W3£ av 



(59) 



Now we relate the average energy of the random motion of particles E &v to 
the distance to Saturn. Here we use results of a hydrodynamic analysis of 
the rings kinetics [29], which from the energy budget of the ring a relation 
between the Kepler frequency Vt = ^GM/r 3 (M - is the mass of Saturn, G 
- gravitational constant, and r - the distance to Saturn), diameter of rings 
particles d, equatorial packing fraction z/ and average square of random 
velocity T : 



fid 



128(1 - e 2 )(l + e)(3 - e) / °° dzu 2 f 3 / 2 
5^ Jo°° dzvf 



(60) 
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where e is the restitution coefficient of the ring particles and v = v(z)jv§ 
and T = T(z)/T are respectively the reduced ^-dependent packing fraction 
and mean square random velocity. Using R av for d and V 2 for T , we recast 
this equation into the form, 

V 2 = Po n 2 Rlu \ (61) 

where the constant po depends on the restitution coefficient and the vertical 
distribution of particles. The same relation as (61), but with a different 
coefficient po, is obtained if one uses a triaxial kinetic approach to the rings 
dynamics [30]. One expects po on the order of 0.1 — 0.001 2 . Using the density 
of Saturn, equal to 700kg/m 3 , the relation R av = 0.506i? c and expressing the 
distance to Saturn as f = r/R s (R s is the Saturn radius), we obtain, 

V = Pl R c u 1 f- 3 / 2 , (62) 

with pi = 2.24 x 10~ 4 v /p^, which together with Eq. 58 yields, 

E agg /E av = p 2 u 2 R- 11 / 3 f 3 , (63) 

where p 2 = 0.0652/p . Exploiting the relations, derived in the previous Sec. 4 
for E agg and E av , expressed in terms of R av , V and r , and Eq. 62 we find, 

E irag /E av = p 3 v 2 r - 5/3 R- 2 f 3 , (64) 

with p 3 = 0.0196/po- Substituting Eqs. 63 and 64 into Eq. (59) we arrive at 
the equation for the upper cutoff radius: 

R c = Bf p exp(af 3 /^) , (65) 

where 

0.583 9/53 24/53 R 36 0.00222 z/g 

B = ^2753 r o ^0 , P = ^ - 0.68 , a = ^~ 

Po b6 Po r ' 

The particular quantities in the last equation are determined by the value of 
Po, which depends on the restitution coefficient of particles and the vertical 
mass and energy distribution in the rings. 

2 For instance, for e = 0.7 and for the Gaussian approximations v(z) = T(z) = 
cxp(— z 2 /h 2 ), where h is the thickness of the rings, the constant po reads, p = 0.0814. 
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Figure 4: The dependence of the upper cutoff radius R c on the distance from 
Saturn r s . Dots represent the observational data taken from table III of [22], 
for the C ring, Cassini Division, and the A ring. Lines - the prediction of the 
theory with a = 12.56. Upper (black) line - the asymptotic theory, Eq. 66, 
the bottom (blue) line - the full theory, that is, the solution of Eq. 65. 



One can find an approximate solution to Eq. 65, depending on the value 
of af 3 . If ct 1//2 f 3//2 ^> 1 (large density limit) the following approximate ex- 
pansion holds 

nA/ 2 f 3 / 2 

R 



yiog(a 1 /2f3/2) _ | log (V/2 f 3/2 _ I loglog ( a l/2f3/2 _ 
a l/2f3/2 



(66) 



where the terms in the "nested" expansion in the above equation are kept, 
until the argument of the "nested" logarithmic function remains positive. 
In the opposite case a 1//2 f 3 / 2 C 1 we obtain (low density limit), 

R c ~ B (l + a/B 2 f 3 " 2/3 + ...)- B f 068 (l + a/B 2 f L64 + . . .) (67) 



In Fig. 4 the dependence of the upper cutoff radius R c on the distance 
from Saturn is shown together with the theoretical predictions of the full 
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theory, Eq. 65, and asymptotic theory, Eq. 66, for the parameter a = 12.56. 
Although the observational data is very scattered, one can still conclude that 
the theory adequately describes the observed behavior of R c . 
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